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Abstract 

A modified version of the interlayer tunneling model, including interlayer single 
particle hopping (ISPH), is considered as a phenomenological model to describe 
cuprate superconductors. The effective ISPH (t±) is taken along with a proba- 
bility factor P, that involves the normal state pseudogap (E g ). This makes t e _^ 
to mimic experimental observations that, ISPH is small in the underdoped regime 
and increases towards overdoping. Within the modified model, we establish the 
absence of bilayer splitting as observed in case of layered cuprates. Transition 
temperature (T c ) and the superconducting gap are calculated. A match, to the 
T-dependent superconducting gap data from experiment, is obtained and high 
values of the ratio of the superconducting gap to T c are recovered. Depending on 
the values of E g , T c as a function of interlayer coupling shows mixed behaviour. 
This is a prediction and can be checked further. 
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Surge of research activity on high temperature superconductors (HTS) in the last few 
years, have clearly established that, the HTS have unusual normal as well as supercon- 
ducting state properties different from those of the conventional BCS superconductors. 
In the superconducting state, the issues related to the superconducting gap, such as the 
symmetry of the gap, ratio of the gap to T c and the temperature variation of the gap, 
are worth addressing. There is growing evidence that the gap of the HTS has d x 2_ y 2 
symmetry with line nodes |E|, §, in contrast to the s-wave gap of the BCS supercon- 
ductors. Similarly, temperature dependence of the superconducting gap and the value 
of the gap ratio, as observed in experiments |3|, |J], are quite different from those of the 
conventional superconductors. 

Endeavour to explain various unusual properties of HTS, resulted in the introduction 
of a number of phenomenological models which are commonly characterized by a BCS 
like gap equation together with well defined quasiparticles in the superconducting state. 
One such model is the interlayer tunneling model originally proposed by P. W. Ander- 
son |y. The interlayer tunneling (ILT) model describes bilayer cuprates where interlayer 
single particle hopping (ISPH) is taken to be inhibited due to correlation effects, but 
tunneling of Cooper pairs between the layers is considered. This pair tunneling amplifies 
the pairing mechanism within a Cu02 layer and an increase in transition temperature 
(T c ) results |J. The ILT model could account for high T c observed in cuprate super- 
conductors, but there are difficulties in recovering very high values of the gap ratio for 
realistic parameters as well as matching the experimental gap variation data within the 
model. At this point, a relevant question to ask is whether the introduction of ISPH 
necessary. As a matter of fact, in band structure calculations, the ISPH matrix elements 
are found to be substantial for bilayer cuprates 0. In addition, experimental observa- 
tions on sufficiently overdoped cuprates find metallic-like behaviour of c-axis resistivity 
and the presence of Drude peak in the c-axis optical conductivity || |J. These features 
immediately imply that, in case of HTS the existence of ISPH cannot be ruled out at 
least in the overdoped regime. 
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Motivated by the experimental observations M, [9j, band structure calculations [|7| and 
the discrepancies reflected in the properties involving the superconducting gap within 
the ILT model, we consider ISPH along with pair hopping to modify the ILT model. For 
underdoped cuprates, the c-axis resistivity shows semiconducting behavior and Drude 
peak is absent in the optical conductivity along the c-axis. These are signatures of 
the absence or marginal presence of ISPH in the underdoped systems and hence, it is 
considered along with a probability factor P such that the effective ISPH is very small 
for underdoping and attains significance for overdoped systems. Pair hopping involves 
the particles which are not taking part in the ISPH. Hence, the probability for pair 



hopping is defined [01], III] as (1 — P) 



By a detailed study of the density of states within the modified ILT model, we address 
the issue of bilayer splitting for layered cuprates. Obtaining the gap equation within the 
mean-field approximation we calculate the superconducting g function of temper- 

ature, and at optimal doping we study the variation of T c with the interlayer coupling. 
Principal results from our calculations are as follows. We demonstrate the absence of 
bilayer splitting within the modified model validating the inclusion of ISPH. Variation of 
T c with interlayer coupling is different for different values of E g . Temperature variation 
of the gap from our calculations has an excellant match to the experimental data 0] and 
ratio of the superconducting gap to T c has high values as in experiments ||. 



The Model : The Hamiltonian of the coupled bilayer complex [10, 11] is written as 



H = E - M)<M - E [V** cWc%c% lC U + h.c 

i,k,ir i,k,k' 

- E [ T f( k ) <$c%c%c§ + h.c] + E [tf(k)cg t cg + /i. C 

i^j,k *^j,k,o- 



This is similar to the one considered by Anderson and coworkers [F2], barring the last 
term accounting for the interlayer single particle hopping. The operator cj^J (cj^) creates 
(annihilates) a fermion in the i-th layer (i=l,2) with momentum k and spin index a. Here 
Vk,k' is the interaction potential forming Cooper pairs in the afr-plane and \x is the chem- 
ical potential. For the afr-plane band dispersion e^, we use the realistic band structure 



obtained from a six parameter tight binding fit [to, ti, £2, £3, £4, £5] = [0.131,-0.149,0.041,- 
0.013,-0.014,0.013] eV to the ARPES data |[J on Bi2212. This six parameter band 
dispersion, used elsewhere |14]| , shows flat bands in the Brillouin zone and the corre- 
sponding density of states (DOS) has a power law singularity known as extended van 
Hove singularity. This is characteristic to the high-T c cuprates Rl5| . The effective ISPH 
matrix element is taken to be t c f(k) = t b ± (k)P, where tj_(k) = t±((cos k x a — cos k y a)/2) 2 
is the k-dependent ISPH adapted from the band structure calculations JF| with t± being 
the bare ISPH matrix element (a is the lattice constant) and P is a probability fac- 
tor determined by the normal state pseudogap found recently in layered cuprates [|16[ . 
Particles that are not taking part in the ISPH, are available for pair tunneling and com- 
plementary probability for each partner of a pair is (1 — P). This fixes the tunneling 
probability of a pair as (1 — P) 2 , and following Anderson || we take the effective pair 
tunneling Tp ff (k) = [(£j_(k)) 2 /|ti|] (1 — P) 2 , where ti is the nearest neighbour hopping 
matrix element of the afe-plane band dispersion (13 . 



Interlayer Single Particle Hopping Probability and the Normal State Pseudogap : In ex- 
periments on layered materials, an important measurable quantity is the c-axis resistivity 
which, when studied as a function of doping (<5), provides usefull information about the 
strength of the effective single particle hopping between the layers. For layered cuprates, 
the c-axis resistivity shows semiconducting behaviour in the underdoped regime JD]] and 



has metallic-like temperature dependence in the heavily overdoped regime || §. This 
has direct bearing on the ISPH and implies that the effective ISPH (t e f) is strong enough 
for overdoping, decreases gradually through optimal doping and becomes weaker towards 
underdoping. The probability factor P, buried in t e f could account for this doping de- 
pendence, when suitably expressed in terms of the temperature and the normal state 
pseudogap E g of cuprate superconductors. The magnitude of E g is found to be van- 
ishingly small for overdoping and it increases towards underdoping ||18|| . Previously, we 
considered an exponential form [[l(], |Tl| for the probability factor, P = exp(—E g /T) and 



found that the experimental gap variation data, in the low temperature region, were 
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not precisely recoverable. We also noticed that a linear T-dependent form for P could 
significantly improve the matching of experimental data. Thus, here we consider a linear 
(low-temperature) T-dependent form (for E g ^> T) of the probability factor as 

T 

P = E g + Te- E s/ T (2) 

where T = with the Boltzmann constant set to unity. Clearly, P — > 1 for E g — > 
(heavily overdoped situation) and P — > as E g — > oo (underdoped situation). Thus P, 
as in Eq(2), consistently mimics the observable doping dependence of the effective ISPH. 
Though a concrete justification of Eq(2) is not possible at present, we present plausible 
arguments in favour of the linear T-dependence of P in connection with the RVB model. 

Within the RVB framework, where spin-charge separation yields spinon and holon 
quasiparticles, Anderson argued that a real hole can come into being in a layer only when 
a holon combines with a spinon and this hole can then hop to another layer. Thus, within 
RVB model, c-axis transport is proportional to the spinon density. Since, in the RVB 
model, the spinon density is proportional to T, hence the linear T-dependence of c-axis 
transport follows. In calculations, we incorporate the momentum dependence of the 
normal state pseudogap, which is found to be of d x 2_ y 2 symmetry [1611, and we take 
E g (k) = E g \ cos k x a — cos k y a\. 

Self Consistent Equations for the Superconducting Gap and the Chemical Potential: 
Within the Hartree-Fock approximation, self consistent equations for the superconduct- 
ing gap Ak and the chemical potential fi (fixing the average number of electrons per site 
n = 1 — 5) are derived for the modified model under consideration. Decoupling of the 
four fermion terms in the Hamiltonian of Eq(l) gives 



i,k,cr 

+ E [t 



(0U0 



+ eff 



(k)cj 



(0t_0") 

kcr C k<T 



i,k 

+ h.c 



k c kj C -k| ~\~ h.C 



(3) 
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where the k-dependent superconducting gap is 

A k = A, k = £ (c% lC V) + Tf (k) (c%c%) (4) 
k' 

which includes pairing in one layer as well as a contribution from pairing in the other 
layer owing to the pair tunneling mechanism. Since the in-plane pairing is identical in 
both the layers, we suppress labelling the gap parameter by layer index. Coherent single 
particle hopping between the layers produces two quasiparticle bands 



K = VK-/i-tf(k)} 2 + A£ (5a) 

and 

Et = V{e k -/i + tf(k)P + A£ (56) 
and the gap parameter turns out to be 

(6) 



l-Tf(k)( X (E k -) + X (E k + ))/2 
where X (E±) = ^t anh (^rV The 

pairing potential can be expanded as Vk,k' = 



2E ± «u. lu , 2 

Vrj^r)^, where rjk is the basis function of the c^ v point group. This makes it possible to 
write A k in terms of a k-independent gap A, the basis function rj^ and the interaction 
parameter V. Finally, the self consistent equation for superconducting gap becomes 

1 i Tfc{x(EZ)+x(E£))/2 

_L = _L^ : '1 (7) 

W N k 1 - Tf(k) ( X (E k ) + X (E+))/2 
and the expression for chemical potential is 

i-^ = i-^E(^-^-^ ff )x(^k)-^E(^-^ + ^)x(^k) (8) 

k k 

Regarding the order parameter symmetry, though the issue is yet to be settled on a 
definite footing, a consensus seems emerging in the many experiments done on layered 
cuprates. The angle resolved photoemission spectroscopy (ARPES) measurements and 
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other phase sensitive experiments indicate the d x 2_ y 2 symmetry of the order parameter 
[0, 0, and we consider the same in our calculations. This means the basis function 
Vk = {cos k x a — cosk y a)/2, which changes sign under tt/2 rotation. The interaction 
parameter V is thus the nearest neighbour and attractive. 

Bilayer Splitting : For bilayer cuprates, the presence of effective single particle hopping 
between the layers should be observable in the single particle density of states (DOS). 
In our calculations, the realistic six parameter dispersion for Bi2212, used to model the 
in-plane band structure, shows a power law singularity in the DOS in the form of a 
single peak. Inclusion of ISPH between the Cu02 layers would split the bands and the 
corresponding DOS will have two peaks. This phenomenon is known as bilayer splitting. 
However, no splitting is observed experimentally in the Bi2212 systems [[HJ] even at 
low temperatures where broadening due to finite lifetime of quasiparticles is expected 
to be small. In Fig.l, we plot the DOS, iV(£) at two different temperatures where 
£k = £k ± £jf(k). Temperature dependence of iV(£) comes only through the probability 
factor P. Parameters for the plots are chosen to be t± = 40 meV and E g = 8 meV, same 
as used to match the gap- variation data in Fig. 3. At low temperatures T = 10 K and 
with no broadening (i.e. broadening parameter T = 0), we find that N(£) has two peaks 
(solid line) separated by an energy ~ 3 meV. However, inclusion of a small broadening 
T = 3meV (dashed line) smears the two peaks and one broad peak appears. At a 
higher temperature T = 50 K, the larger peak separation energy ~ 21meV is due to 
the enhanced probability of the interlayer single particle hopping. Inclusion of nonzero 
T (dashed lines) broadens the peaks and the value of T needed for complete smearing of 
two peaks is T > lOmeV. This is demonstrated in the lower panel (T = 50 K) of Fig.l. 

It should be mentioned here that, the temperature dependence of the broadening 
parameter for layered cuprates is noted to follow || the relation V (meV) ~ To + 
2.5T C (T/T C ) 3 , where r is the zero temperature value. Although high values of T ~ 
lhmeV is observed for Bi-cuprates ||, in our case a T ~ 7 — 8meV could be enough 
to remove the splitting at all temperatures. Moreover, at low temperatures (T = 10 K) 
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the peak separation is small ~ 3 meV (solid line for T = 10 K in Fig.l), whereas in the 
ARPES experiments the energy resolution is far beyond this level (resolution function 



FWHM~ l9meV) |20| and the small peak separation would probably remain unre- 
solved. Thus, in the cleanest sample and at very low temperature, even if one ideally 
takes To — * 0, the bilayer splitting in Bi-cuprates would remain unobservable by exper- 
iments. Finally, it has been noted that YBCO systems have a smaller Tq ~ 1.5 meV 
[H] and also have orthorhombic distortions, which might explain why splitting could 



possibly be observed in YBCO plf in contrast to the Bi-cuprates. 



So far, we have discussed the issue of bilayer splitting by calculating the normal state 
density of states. But, in actual ARPES experiments, bilayer splitting is addressed by 
studying the photoemission intensity curves in the superconducting state. Within the 



modified ILT model, we have calculated the photoemission intensity curves [11 1 in the 



superconducting state and found that, with an energy resolution much better than that 



in actual experiments |20fl and with a very small boradening T ~ 1 meV, bilayer splitting 
remains unobservable for temperatures ranging upto 50 K. In other words, the inclusion 
of an effective ISPH in the modified ILT model is in perfect tune with the experimental 
observations. 

Issues related to the Superconducting Gap : The superconducting gap and other quanti- 
ties are calculated by simultaneous solutions of the self consistent Eq(7) and Eq(8) for 
the gap and the chemical potential respectively, implementing numerical techniques. 

In our calculations, doping is kept fixed at the optimum level (S opt ) at which T c 
attains its maximum value T c m . Position of 5 opt is found to have a slow dependence 
on the value of E g considered. A study of T c versus 5, for different E g , shows that 
with Eg < 5 meV, T c has a two-peak structure, whereas for higher E g a one-peak bell- 
shaped form is recovered consistent with that in layered cuprates. Thus, we enforce the 
condition E g > 5 meV. Value of bare single particle hopping matrix element is taken to 
be t± = 40meV in accordance with the band structure calculations suggesting ~ 
0.2 — 0.3. Interaction parameter is taken to be V = 80meV which makes T c m ~ 100 K 



8 



for Eg ~ T™. This value of V is kept fixed throughout in this communication. 

In Fig. 2 (a) we plot T c m as a function of t±. Different curves, labelled by capital 
alphabets, correspond to various E g values as, A:6, B:10, C:15 and D:20 in meV. Clearly, 
variation of T c m with t± may take different shapes depending on the values of E g . For 
small Eg, the suppression of coherent single particle hopping between the layers is weak 
and T c m decreases with t±, whereas large E g makes ISPH less probable and in turn the 
pair tunneling probability gets enhanced, which results a rise in T c m with t±. In the 
Anderson limit (E g = oo), T c m grows rapidly with t±. 

Application of pressure on a sample along c-axis could reduce the out-of-plane lattice 
constant which results in an increase of the interlayer coupling. In the plot of T c m versus 
t± in Fig. 2 (a), we find that, for large E g (underdoped materials), T c m increases with t±, 
for small E g (overdoped materials), T c m decreases with t± and at some intermediate E g , 
rpm rerri ains unchanged. Thus, within our model, c-axis applied pressure would result, an 
increase of T c for underdoping, a decrease of T c for overdoping and an unaffected T c for 
some intermediate doping. Similar varied behaviours of T c are seen for differently doped 



cuprates under hydrostatic or uniaxial pressure ||22|| . However, it may be mentioned that 
hydrostatic pressure could also inflict changes in the in plane parameters. Hence, we 
leave this result of various T c m variation with t± as a prediction which could be checked 
by further experiments. 

Variation of the zero temperature superconducting gap (A^ ax (0)) as a function of 
t± is shown in Fig.2(b). The E g values corresponding to different curves are same as 
used in Fig. 2(a). Here, irrespective of the values of E g , A^ ax (0) increases with t± in the 
same fashion for every curve. 

A plot of T c m as a function of E g , for different t±, is given in Fig. 3. As is obvious, 
T c m increases with increasing E g due to the progressive enhancement of effective pair 
tunneling probability, caused by E g . As is seen, at some E g close to ~ 12meV, T™ for 
all the curves are same, implying that the transition temperature is independent of t±. 
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In fact, we crosscheck that, T c m as a function of t±, for E g m 12mev, is nearly flat. 

Next, we come to the issue of the temperature dependence of the ci-wave super- 
conducting gap. A plot of scaled supercondcting gap A^ ax (T) / A^ ax (0) as a function 
of reduced temperature T/T c is given in Fig.4. Solid line is the BCS-form of gap- 
variation and solid squares are experimental data on Bi-cuprate [|J. Except at the edges 
(T/T c = 0&1) the locus of the experimental data is well below the BCS-curve and 
the absence of low temperature saturation in the experimental data stands out as an 
anomalous characteristic of high-T c cuprates. The dashed lines are from our calculations, 
but, the one very close to the BCS-curve represents the gap variation for the original 
ILT model (E g = oo). The dashed curve with E g = 8meV is within the modified ILT 
model which matches experimental data very well. To the best of our knowledge, this 
is the first ever calculations where precise matching of experimental gap-variation data 
is obtained. An important and related quantity is the ratio of the zero temperature 
superconducting gap to T c . In experiments on layered high-T c cuprates, high values of 
2A^ ax (0)/T c ~ 5 — 7 is observed |§, compared to the BCS superconductors where it 
is ~ 3.5. In the inset of Fig.4 we plot the gap-ratio as a function of t± for the same 
parameters used to fit the experimental gap- variation data. The gap-ratio increases with 
t± and within the realistic range of t± (~ 30 — 45 meV) we recover the ratio ~ 5 — 7. To 
be precise, at t± = 40meV, we get 2A"^ ax (0)/T c = 5.8. Note that, for same parameters 
and with E g = oo (Anderson limit), we find 2A^ ax (0)/T c = 4.4. High value of the gap- 
ratio and its increasing trend with t± within the modified model, could be understood 
as follows. 

For any finite E g , when T — > 0, the probability factor P —>■ 0, ISPH is completely 
blocked and A^ ax (0) grows with t± since t± acts only to enhance T£ s . But, for any finite 
T > the probability factor P / 0, blocking of the effective ISPH is not complete and 
one gets a transition temperature T C (P ^ 0) < T C (P = 0) resulting in a high value of the 
gap-ratio. Regarding the increasing trend of gap-ratio, consider the situation when T c m 
decreases or saturates with t± for small E g (as in Fig. 2 (a)). But notice that, A^ ax (0) 
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always rises with t± (see Fig.2(b)) yielding the observed result of the gap-ratio. Even 
for the case when T c increases with t± (for high E g ), note that its increase is slower 
than A^ ax (0) which effectively increases the gap-ratio with t±. For a fixed value of t±, 
one gets an increasing gap-ratio with decreasing E g which could also follow from similar 
arguments. 

To conclude, we have considered a modified interlayer pair tunneling model where 
pair tunneling is accompanied by interlayer single particle hopping. Within the model, 
absence of bilayer splitting is established corresponding to experimental situation. Dif- 
ferent types of variation of T c as a function of interlayer coupling, for different values 
of the pseudogap, comes out as a prediction that corresponds to T c variation under 
c-axis pressure. Temperature dependent superconducting gap from our calculations cor- 
rectly reproduces the experimental data, and for realistic parameters, high values of the 
gap-ratio are obtained. 
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Figure captions: 



Fig.l. Single particle density of states N(£) is plotted at two different temperatures 
T — 10 K and T = 50 K. Values of broadening parameter V used are given in the 
figures. The full bandwidth of energy £ is not shown, since the peaks in the DOS are 
only of concern for the purpose of addressing the issue of bilayer splitting. 

Fig. 2. (a) Mean-field transition temperature (T c m ) at optimal doping as a function of 
bare interlayer coupling t± for different E g represented by alphabetic labels as A:6, B:10, 
C:15, D:20 in meV. The interaction parameter is V = 80meV. 

(b) A plot of the maximum value of the zero temperature superconducting gap. Values 
of Eg for different curves are as shown and the alphabetic labels corresponds to same 
Eg values as in (a). 

Fig. 3. Transition temperature as a function of the normal state pseudogap magnitude 
Eg for three different values of t±, as shown in the figure. Doping is set to the optimum 
level (for E g = 8 meV) and V = 80 meV. 

Fig. 4. Temperature variation of the scaled superconducting gap (A^ iax (T)/A^ iaa: (0)) as 
a function of the reduced temperature (T/T c ) is presented. Solid line corresponds to the 
conventional BCS superconductors and solid squares are experimental data from Ref. @]. 
Dashed lines are from our calculations for different E g shown. [Inset: A plot of the ratio 
A^ ax (0)/T c as a function of t±. The parameters are V = 80meV and E g = 8meV that 
yielded match to the experimental gap-variation data as in the main figure.] 
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